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X-RAY MONITORING 

The present invention relates to the monitoring of mixtures, and in 
particular to real time monitoring using X-ray scanning. It has application 

5 in a number of fields were the monitoring and imaging of mixtures is 
required, and is particularly applicable to the measurement of the flow 
rates of the different fractions in moving mixtures of fluids, such as the 
measurement of oil flow from oil wells. It is also applicable to slurry 
flow monitoring and measturement, the monitoring of fluidization 

10 processes, and in the monitoring of mixing processes. 

Typically, an oil well is one of a number of wells that collectively form 
an oil field. As well as extracting oil. water and gas from an individual 
oil well, it is sometimes necessary to recycle water down an oil well back 
15 into the ground. Typically, this pumped water is used to facilitate 
diffusion of oil through the porous rocks in the grotmd towards a well (or 
wells) for subsequent collection. By careful design of pumping and 
extraction rates* it is possible to maximise the yield of oil from the field. 

20 Often, the output from several extracting oil wells are joined into a single 
pipeline for subsequent downstream processing. To optimise the 
production process, it is necessary to know the fraction of oil and water 
in the liquid phase and also to know the volumes of oil and water 
produced by each well in the field. By combining this information with an 

25 understanding of the geophysics and seismology data of the field itself, it 
is believed to be possible to improve on the quantity of oil produced from 
the field and to reduce the cost of production. Such measurements of oil 
and water phase fraction and velocity require instrumentation to be placed 
at the well head. 

30 
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Known instrumentation for this purpose relies on measurement of linear 
attenuation coefficient of the flow using a dual-energy gamma 
measurement once the flow has been homogenised by use of a mechanical 
system. A second known approach is to measure electrical properties of 
5 the fluid (including permittivity and conductivity) and from this to infer 
phase fraction of oil and water. 

The present invention provides apparatus for monitoring in real time the 
movement of a plurality of substances in a mixture, the apparatus 
10 comprising an X-ray scanner arranged to make a plurality of scans of the 
mixture over a monitoring period to produce a plurality of scan data sets, 
and control means arranged to analyse the data sets to identify volumes of 
each of the substances and to measture their movement. 

15 Preferably the apparatus is arranged, on each scan, to produce a data set 
relating to a layer of the mixture. More preferably the control means is 
arranged to define a plurality of volume elements in said layer and to use 
a measure of the X-ray attenuation in each of said volume elements to 
form the data set. Still more preferably the control means is arranged to 

20 use the data sets to determine the amount of at least one of the substances 
in said layer. 

Preferably the control means is arranged to use the data sets from each of 
the scans to determine a time averaged value of the amount of said at 
25 least one substance. 

Preferably the scanner is arranged to produce data sets relating to a 
pltirality of layers of the mixture, the layers being in different positions 
from each other. More preferably the control means is arranged to use the 
30 data sets relating to said plurality of layers to measure movement of at 
least one of the substances. Still more preferably the control means is 
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arranged to track the movement of regions of said substance through the 
plurality of layers to determine a flow velocity of said substance. 

The control means may be arranged to measure the movement of a region 
5 of a first one of the substances, to determine a measure of the buoyancy 
of said region relative to at least one other substance, and to measure the 
movement of said at least one other substance using the movement of said 
region and said buoyancy. 

10 Preferably the control means is arranged to define a model for calculating 
a parameter of movement of the substances on the basis of a number of 
variables, to produce a measured value of the parameter from the scan 
data sets, and to determine at least one of said variables from the 
measured value and the model. 

15 

The control means may be arranged to determine a flow rate of at least 
one of the substances, the flow rate being defined as the amount of said 
substance flowing through a predetermined region in a predetermined 
time. 

20 

Preferably the control means is arranged to analyse a scan data set in two 
stages, one stage providing a lower spatial resolution and higher contrast 
resolution than the other. In this case the control means is preferably 
arranged to use the higher spatial resolution analysis to identify volumes 

25 of a first of said substances and to use the higher contrast analysis to 
distinguish between volumes of two further substances. For example the 
control means may be arranged to use the high spatial resolution analysis 
to adjust a measure of X-ray attenuation, of volume elements defined in 
the low spatial resolution analysis, to account for the presence in said 

30 volume elements of the first substance. This is particularly useful where 
the mixture contains two substances of similar X-ray attenuation 



wo 2004/097344 



PCT/GB2004/001751 



4 

coefficient, for example both being in a liquid phase or both being in a 
solid phase, such as oil and water, and another substance of a 
significantly different X-ray attenuation coefficient, for example being of 
a different phase e.g. a gas, such as air. 

5 

The scanner may be arranged to be placed around a pipe or other conduit 
to measure the movement of the substances through the pipe or conduit. 
Altemativcjly the scanner may be placed adjacent to or around a container 
in which the mixture is undergoing a fluidization or mixing process. 

10 

The apparatus may further comprise display means arranged to display an 
image of the mixture controlled by the control means. For example the 
display means may be arranged to display a video image of the mixture. 

15 The present invention further provides a method of monitoring in real 
time the movement of a plurality of substances in a mixture, the method 
comprising making a plurality of X-ray scans of the mixture over a 
monitoring period to produce a plurality of scan data sets, and analysing 
the data sets to identify volumes of each of the substances and to measure 

20 their movement. 

Preferred embodiments of the present invention will now be described by 
way of example only with reference to the accompanying drawings in 
which: 

25 

Figure 1 is a transverse section through an X-ray scanning 
apparatus according to a first embodiment of the invention; 



30 



Figure 2 is a longitudinal section through the apparatus of Figure 
1; 
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Figure 3 shows part of a map of a volume produced by the 
apparatus of Figure 1; 



Figure 4 is a graph showing how the error in the volume 
5 measurement using the apparatus of Figure 1 decreases with time; 

Figure 5 shows how flow velocity is determined using the 
apparatus of Figure 1; and 

10 Figure 6 is a schematic section through a scanning apparatus 

according to a second embodiment of the invention. 



Referring to Figures 1 and 2 an X-ray scanner 8 according to a first 
embodiment of the invention comprises an annular multi-focus X-ray tube 

15 10 and an annular segmented X-ray sensor array 12 located radially 
inwards of the X-ray tube 10. The tube 10 comprises a number of X-ray 
sources 14 spaced around the tube 10, each X-ray source 14 comprising a 
target 16 at which electrons are directed from an electron source 18, The 
sensor array 12 is cylindrical in form and comprises a number of sensor 

20 elements 20 arranged in a nxxmber of adjacent rings 22a, 22b. The X-ray 
sources 14 are just above the top of the sensor array 12 and each source 
14 is arranged to direct X-rays over the part of the sensor array 12 
nearest to it, and towards a wide area of the sensor array 12 on the 
opposite side of the scanner axis Z. The scanner 8 is placed around an oil 

25 pipe 24 so that the pipe 24 lies on the scanner axis Z. The X-ray system 
would normally be operated at 100 - 150 kVp tube voltage and 10 - 50 
mA beam current. 

A control unit 26 is connected to each of the electron sources 18 so that it 
30 can control which X-ray source 14 is active at any time, and also to each 
of the sensor elements 20 in the sensor array 12 so that it can determine 
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the strength of the X-ray radiation reaching each sensor element 20 at any 
time. The control imit 26 controls the X-ray sources to scan repeatedly at 
a. scan period of between 0.1 and 10ms per scan to acquire scan data, 
which can be analysed to provide quantitative outputs, or to provide 
5 images of the mixture on a screen 27. 

During each scan, each of the X-ray sources 14 emits once in turn, and 
for each source an image data set is formed for each of a number of 
image rings, one image ring for each ring 22a, 22b of the sensor array 

10 12. Each image ring is at a slight angle, in this case between 1 and 5 
degrees from the normal to the pipe axis Z, due to the axial offset 
between the sources 14 and the sensor array 12. The control unit 26 
processes the image data sets from each of the image rings produced in 
the scan to produce a tomographic image data set made up of a number of 

15 plane image data sets for respective image planes 28a, 28b spaced along a 
length of the pipe 24 as shown in Figure 5. Typically there are between 2 
and 14 planes at a plane spacing of 1 to 2mm. Each plane image data set 
includes a grey level for each pixel of the image, which represents the X- 
ray attenuation coefficient of each of a number of corresponding voltune 

20 elements (voxels) in the respective 'plane' or layer of the imaged volume. 
The layer has a thickness equal to the thickness of one volume element, 
making the plane image essentially two dimensional. 

The tomographic image data sets are then processed in two stages, the 
25 first to determine the relative volumes of the various fractions in the pipe, 
in this case oil, water and gas, and the second to determine the absolute 
velocities of the oil, water and gas phases. The volume flow rate for each 
phase can then be determined as the product of phase fraction with phase 
velocity. 



30 
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To determine the phase fractions two image reconstructions are generated 
for each image plane 28a, one at low spatial resolution but good contrast 
resolution and one at high spatial resolution but lower contrast resolution. 
Both image reconstructions use the same initial projection data set. The 

5 area of each image pixel corresponds to the cross sectional area of the 
corresponding voxel, and both reconstructions have the same slice 
thickness. Figure 3 shows an area of the image plane 28a divided at both 
high and low resolution. Typically, the high-resolution image will be 
reconstructed with smaller pixels 30 of 2 nmi pixel dimensions, while a 

10 low-resolution image may be reconstructed with larger pixels 32 of up to 
10 mm pixel dimensions. In the example shown the smaller pixels 30 are 
of 2mm dimension and the larger pixels 32 are of 6mm dimension. In 
each image, each pixel will have a grey level indicative of the attenuation 
coefficient of the substance or substances in the corresponding volume 

15 element. For any volume elements which contain only one fraction, the 
grey level will be at one of three possible levels corresponding to oil, 
water and air. For any volume element containing more than one fraction, 
the grey level will be at a level between these three levels. 

20 Firstly, the high resolution image is segmented, using a suitable binary 
segmentation process, to determine which of the smaller pixels 30 
represent the gas phase, and which represent the liquid phase. This relies 
on the fact that the gas phase reconstructs back to a very different signal 
level from that of the liquid phase, so such binary segmentation is 

25 straightforward. The number of pixels representing the gas phase 
therefore gives a measure of the gas fraction in the image layer. In Figure 
3, the area 34 is representative of a gas bubble which fills one of the 
large pixels 32 and a further 9 of the small pixels 30, spread between 3 of 
the large pixels. 
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Next, any gas voxels that are identified from the high spatial resolution 
scan are used to apply a partial volume correction to the low spatial 
resolution scan as shown in Figure 3. This correction identifies that a 
fraction of the low spatial resolution pixel is filled with gas, and that 

5 therefore the measured grey level is lower than possible for a liquid phase 
only voxel. For example the middle large pixel 32a includes 3 out of 9 
small pixels 30 representing gas voxels. These are the shaded pixels 30a« 
A linear volume correction is applied to correct the grey level of the large 
pixel 32a to accotmt for the gas partial fraction within the corresponding 

10 voxel: 

true = — actual 

where true = partial volume corrected grey level, Vgas = gas fraction of 
large pixel and actual = actual reconstructed grey level. 

15 When these corrections have been made, the large pixels 32 are 
segmented on the basis of the corrected grey levels into those 
representing voxels in which gas predominates, those representing voxels 
in which water predominates, and those representing voxels in which oil 
predominates. This therefore provides a measure of the instantaneous 

20 volumetric fractions of air, oil and water in one image plane 28a. It also 
provides a tomographic image of the mixture flowing through the pipe 24 
which can be displayed on the screen 27. 

Since the flow is continuously changing, it is necessary to repeat this 
25 measurement multiple times to achieve a time-averaged measurement of 
the volumetric fractions of the three phases. This is shown graphically in 
Figure 4 which shows a number of consecutive instantaneous phase 
volume fraction measurements 40 which are spread around the mean 
fraction. As the number of measurements increases the uncertainty in the 
30 mean volume fraction decreases. Measurement uncertainty depends on the 



wo 2004/097344 



PCT/GB2004/001751 



9 

number of voxels per image and the photon statistics within a voxel. 
Typically, uncertainty in phase fraction within a single scan is at the 5% 
level. Uncertainty reduces to less than 1% after averaging of information 
from multiple data sets. 

5 

Also as a new image is produced for each scan, and the scanning period is 
typically between 0.1 and 10 ms as indicated above, the images produced 
can be combined to provide real time tomographic imaging of the flow in 
the pipe. 

10 

Referring to Figure 5, flow velocity can be determined by applying cross- 
correlation methods between identified features in the stack of image 
planes 28a, 28b that are generated along the length of the pipe 24 at 
different times. This allows the speed of movement of those features 

15 along the pipe to be determined. In the example shown, the gas bubble 34 
shown in Figure 3 moves up the pipe 24. At a reference time t=0 the 
bubble 34 reaches one of the image planes 28d. At the time of the next 
scan, at time t=:0 + 5 the bubble 34 is partly in the image plane 28d, and 
partly in the adjacent image plane 28c. At the next scanning time t 

20 =0 + 28 the bubble 34 is also partly in the next image plane 28b. Using 
cross correlation between the images in the planes 28a, 28b, 28c, 28d 
allows the bubble 34 to be identified as the same feature, and its position 
at the time of each scan to be determined. Using the change of position 
between scans and the scan frequency, the velocity of the bubble 34 along 

25 the pipe can be determined. The velocity of movement of the air, together 
with the volumetric fraction of the air obtained as described above, allow 
determination of the volumetric flow rate of air through the pipe 24, that 
is the volume of air flowing through the pipe per unit time. 

30 In some simple flows it will be possible to use the technique described 
above to identify features of each of the oil, water and gas phases in the 
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pipe and to measure their individual velocities. In other types of flow 
different techniques are required to determine the flow rates of all of the 
phases. 

5 One technique is to include in the calculation the buoyancy of the various 
fractions. In the example shown in Figure 5, if the gas bubble 34 is 
surrounded entirely by oil, then its speed of motion will be determined 
partly by its buoyancy and partly by the velocity and viscosity of the oil. 
As the approximate viscosity of the oil and the approximate densities of 

10 oil and water are known, measuring the speed of the bubble 34 can be 
used to measure the velocity of the oil. 

Alternatively the density of each of the phases can be measured. The 
density of each phase is very closely related to its electron density, which 
15 determines the X-ray attenuation as measured by the scanner. Therefore 
the measured attenuation of each voxel can be used to determine the 
density of the fluid in it. The measured densities can then be used to 
determine the buoyancy of each of the fractions. 

20 The viscosity of each of the phases can also be measured by measuring 
the rate of change of shape of the features of that phase. For substances 
of high viscosity, such as oil, the rate of change of shape will be 
relatively slow, whereas for substances of lower viscosity, such as water, 
the rate of change of shape will be relatively fast. 

25 

For complex flow patterns mathematical modelling of the flow can be 
used to determine the speeds of flow and flow rates of the different 
phases. A model can be built up which will calculate the size, shape and 
velocity of the regions of each fraction based on the flow rates of each of 
30 the fractions and other variables such as the pipe size, pipe orientation, 
and temperature. The known variables such as pipe size and temperature 
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are input to the model. The size of any regions of each phase that can be 
measured, and the flow velocities of any of the phases that can be 
measured are matched to find a best fit with those calculated using the 
model for different flow rates using a *least squares' or similar method to 
5 determine the actual flow rates. By repeating the best fit process the 
results produced can be averaged to improve the accuracy of the results. 
This modelling technique is particularly suitable for monitoring turbulent 
flow, or flow with high gas fractions where reverse flow can occur. 

10 Referring to Figure 6 in a second embodiment of the invention a bed 100 
of powdered substances is fluidised by the passage of a gas 102 through 
it. The gas 102 acts as a catalyst for a chemical reaction. A scanner 104 
operating in a similar manner to that of the first embodiment is placed 
arotmd the fluidised bed 100 to produce real time tomographic imaging of 

15 the fluidised powder. This can be used to monitor and analyse the 
fltudisation to ensure that it is functioning as required, using a control 
unit 106. 

It will be appreciated that the method described above can be used for 
20 real time imaging, or measuring the movement of, a variety of substances 
in a variety of conditions. For example as well as monitoring or 
measuring liquid and gas flow, it can be used to measure powder flow, in 
which powder is carried in a gas, and slurry flow in which a powder is 
carried in a liquid. It can also be used to monitor other processes such as 
25 fluidisation in which gas is passed through a powder or granular solid, for 
example to. encourage a chemical reaction, and mixing processes to 
monitor how well different substances are being mixed together. 
Depending on the substances being monitored various other parameters 
can be analysed from the data, such as the size and velocity of gas 
30 bubbles, of droplets or regions of liquid, or of solid particles, or the 
degree of mixing of the substances. 



